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Abstract 



The possibility of introducing non-perturbative infrared effects leading to a modification of the long 
distance behavior of gauge theories through a soft breaking of the BRST invariance is investigated. The 
method reproduces the Gribov-Zwanziger action describing the restriction of the domain of integration 
in the Feynman path integral to the Gribov region and a model for the dynamical quark mass generation 
is presented. The soft symmetry breaking relies on the introduction of BRST doublets and massive 
physical parameters, which allow one to distinguish the infrared region from the ultraviolet one, within 
the same theory. 
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1 Introduction 



The task of introducing non-perturbative effects is one of the most difficult challenges in quantum field 
theory. In this letter we propose to introduce non-perturbative infrared effects in a way that preserves 
the basic properties of continuum field theory, namely: locality and renormalizability. The idea relies on 
the possibility of introducing in a controllable way local terms which give rise to a soft breaking of the 
BRST symmetry, meaning that the dimension of these breaking terms is smaller than the space-time 
dimension. Such soft terms can be neglected in the deep ultraviolet region, where one recovers the notion 
of exact BRST invariance. However, in the low energy infrared region, the soft breaking terms can- 
not be neglected, and produce a quantitative modification of the large distance behavior of the theory. 
Furthermore, these terms are introduced in a way that preserves the nilpotency of the BRST operator 
and doesn't modify its cohomology. As a consequence, the space of the local observables of the theory, 
identified with the cohomology classes of the BRST operator in the space of local field polynomials, 
remains unchanged. In other words, assuming that the relevant correlation functions of the theory are 
of the type (Oi(xi)....0„(x„)) with Oi(xi), 0„(x„) field polynomials belonging to the cohomology of 
the BRST operator, we attempt at introducing nonperturbative effects which might modify the infrared 
behavior of (Oi(xi)....0„(a;„)), while leaving unchanged its behavior in the deep ultraviolet region. Thus, 
the same correlation function can display a different behavior, according to the region which is being 
considered. Our aim here is precisely that of keeping the same set of observables of the starting theory, 
while modifying the large distance behavior of their correlation functions. 

The first step in order to achieve this goal is that of introducing a suitable set of additional fields, 
assembled as BRST doublets. Such fields do not modify the cohomology of the BRST operator. Thus, 
their introduction leaves unchanged the original content of observables of the theory. This follows from 
the fact that BRST doublets always give rise to invariant terms which can be cast in the form of exact 
BRST variations, thus yielding vanishing cohomology [1]. 

Secondly, these BRST doublets are coupled to a set of dimensionfull parameters in such a way that 
the initial BRST invariance of the theory can only be softly broken by terms which can be kept under 
control at the quantum level. This nontrivial requirement can be fulfilled thanks to the BRST doublet 
structure of the fields coupled to these dimensionful parameters, which determine suitable linearly broken 
Ward identities, as well as generalized Slavnov-Taylor identities, ensuring the renormalizability and the 
stability of the theory in the presence of the breaking. 

The last step is that of requiring that these parameters are not treated as free parameters of the theory. 
Instead, they are determined in a self-consistent way at the quantum level by imposing appropriate gap 
equations, refiecting the possible non-perturbative condensation of the local operators coupled to them. 

This allows one to express them in terms of the original coupling constant g of the theory, exhibiting 
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the typical non-perturbative behavior ~ e ^ . Albeit we cannot specify in full generality the explicit 
form of these gap equations, it is worth underlining that the dimension of the operators coupled to these 
parameters is smaller than the space-time dimensions. The analytic study of the possible condensation 
of these operators is simplified as compared to the case of an operator of maximum dimension. In fact, 
several techniques are nowadays at our disposal in order to face the study of the condensation of local 
operators of lower dimension: effective potential, local composite operator technique, the renormalization 
group equations, variational principles as the minimal sensitivity principle, the fastest apparent conver- 
gence criterion, etc. Let us mention, for example, the case of the condensation of the dimension two 
mass operator ((^^^ — acc)) in the Yang-Mills theory expressed in the maximal Abelian gauge. The 
condensation of this operator, see ^ and references therein, provides a non-perturbative dynamical mass 
for off-diagonal gluons, a feature of great relevance for the Abelian dominance hypothesis for the dual 
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superconductivity picture of color confinement. 



One can certainly question the meaning of introducing a soft breaking of the BRST symmetry, since the 
latter is known to be related to physical properties of the theory, such as the unitarity and the characteri- 
zation of the physical subspace. The justification is for models displaying non-perturbative phases, which 
require a better knowledge of the infrared sector. This is the case, for example, of the QCD confining 
non-abelian gauge theories. As is well known, the particle interpretation of QCD in terms of partons - 
gluons and quarks- is lost in the low energy region. Unitarity cannot and shouldn't be defined for the 
gluon and quark sector in the infrared region, where confinement sets in, and gluons and quarks do not 
correspond anymore to the physical excitations of the theory. Rather, the physical states are given by 
colorless bound states like baryons, mesons and glueballs. Moreover, the framework that we introduce 
here can also include the description of non-perturbative long distance effects in supersymmetric as well 
as in topological theories, both by taking into account the effect of the Gribov phenomenon and a soft 
breaking of the super symmetry, as in the Seiberg-Witten theory [3]. 

As a last remark, we would like to underline that the soft breaking of the BRST symmetry intro- 
duced here is meant to be an explicit breaking, as opposed to the possibility of achieving a spontaneous 
symmetry breaking. Indeed, the concept of spontaneous BRST symmetry breaking would imply the 
existence of Goldstone massless modes, which would give rise to a rather different framework. 

The paper is organized as follows. In Sect. 2 we give an overview of the introduction of the BRST 
doublet fields and of the related dimensionfull parameters. We point out that the presence of the soft 
breaking of the BRST invariance ensures that the aforementioned dimensionfull parameters are physical 
parameters of the theory and not unphysical ones, as it would be the case of gauge parameters. An 
elementary algebraic proof of this property is given. Thus, in the infrared domain, the values of the 
correlation functions of the theory will explicitly depend on these parameters, which would themselves be 
related to A.qcd by generalized gap equations. In Sect. 3, we show that the method yields the so-called 
Gribov- Zwanziger action which implements the restriction of the domain of integration in the path in- 
tegral to the Gribov region, plus a relevant soft term. In Sect. 4 we present a model for the dynamical 
quark mass generation. 

2 Soft breaking of the BRST invariance 

We start by considering a quantized action 5inv(0), endowed with a set of fields generically denoted by 
(j), with the following properties: 

• the action exhibits exact BRST invariance 



• 5inv is multiplicatively renormalizable leading, in particular, to a consistent perturbation theory in 
the ultraviolet region. 

A typical example of such an action is provided by non-abelian SU {N) Euclidean Yang-Mills theories, 
namely 



s5'inv((/') = 



(1) 




= 




inv — 



SyM + Sgf 



(3) 
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where Sym stands for the Yang-Mills action 

Sym = \ I d'^xF^.F^,, (4) 

and Sgf is the gauge-fixing term of the Landau gauge 

5',/ = J d^x (i6"a^^;i + c'^a^Z)^^c^) , (5) 

with 6° being the Lagrange multiplier field enforcing the Landau condition, = 0, and (c°',c^) the 

Faddeev-popov ghosts. As required, expression ^ is renormalizable, while displaying BRST invari- 
ance 

We will now extend this local action for tentatively introducing non-perturbative infrared effects within 
the context of local quantum field theory, which will lead to a modification of the long distance behavior 
of the theory, while leaving unaffected its ultraviolet behavior. To do so, we introduce a set of fields, 
generically denoted as (a, /3) , which transform as BRST doublets 

sa = P , 

sP = . (6) 

Such fields give rise to a trivial cohomology, meaning that any local BRST invariant functional F{a, (3) 
has necessarily the form of an exact BRST term 

sF{a,P) = ^ F{a,p) = sF{a,P) , (7) 

for some local F{a,(3). 

After the introduction of the fields (a, /?) and of a corresponding BRST invariant term describing their 
propagation, 

S^l{a,P) = sS{a,§) (8) 

one defines a set of dimensionfull parameters (a) and a local term So- (a, a, /?, (p) which gives rise to a soft 
breaking of the BRST invariance, as expressed by 

s (5inv(0) + S^^^{a,P) + S^{a,a,P, </-)) = sS^ = aA{a,P, . (9) 

Since (a) are dimensionfull parameters, the breaking term A(q, /?, (f>) is an integrated local field polyno- 
mial whose dimension is smaller that the space-time dimension, i.e. it is a soft breaking. As such, it 
can be neglected in the deep ultraviolet region. Moreover, taking into account that the additional fields 
(a, (3) form BRST doublets, it follows that the ultraviolet behavior of the model is left unmodified. 

The introduction of the soft parameters (ct) has to be done in a way that preserves the multiplicative 
renormalizability of the modified theory. This physical requirement allows one to carry out calculations 
in a consistent way. Furthermore, it puts severe constraints on the way the soft parameters are intro- 
duced. The simplest and consistent way is to couple the additional doublet fields (a, f3) linearly to the 
original fields cj), so that the resulting breaking term is a quadratic term in the fields. As a consequence, 
the propagator of the i;^-field gets modified in the infrared by the soft parameters entering the quadratic 
coupling. In turn, this will affect the infrared behavior of the correlation functions of the theory, as can 
be seen, in particular, from the modifications at small momentum that are brought to the propagators of 
partons. Coupling the fields in that way has important consequences on the Ward identities of the theory. 
In fact, a quadratic term means that the additional fields (a, /5) couple linearly to the original fields (p. 
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As we shall see in the examples discussed in the next sections, this feature enables one to write down a 
suitable set of linearly broken Ward identities, controlling the renormalizability of the theory. In sum- 
mary, the introduction of the infrared soft breaking can be done in a way in which the resulting action is 
stable against quantum corrections. This will ensure that no additional parameters have to be introduced. 

The presence of the breaking term A(a, /3,(/>) enables one to establish that the dimensionfull param- 
eters (ct) are physical parameters of the theory, as for instance the Gribov parameter 7 [H |5] of the 
Gribov-Zwanziger action. In fact, taking the derivative of both sides of eq.® with respect to o_, one gets 

.^.A. (10) 
from which it is apparent that cannot be cast in the form of an exact BRST variation, namely 



da 



^sA, (11) 



for some local A. Equation (jlip ensures that the parameters (a) are physical parameters of the theory. 
As such, they will influence the expressions of the correlation functions of the operators belonging to the 
cohomology os s, modifying their long distance behavior. We see therefore that the soft breaking term 
A plays a key role in order to introduce nontrivial parameters. Suppose in fact that, instead of giving 
rise to a soft breaking, the term Sa{a,a, l3,(j)) would be left invariant by the BRST operator, i.e. 

sSa{a,a,^,(f)) = . (12) 
Therefore, owing to the doublet structure of (a, /3), a local functional should exist such that 

Sf^{a,a,f3_,(l)) = sSa , (13) 

from which it would follow that 



da da ' 



(14) 



which would imply that (a) would be unphysical parameters. As a consequence, the correlation functions 
of the observables of the theory would be independent from (a), meaning that the introduction of these 
parameters would have no effects on the theory. 



3 The example of the Gribov-Zwanziger Lagrangian 

As a first example, we give a derivation of a refined version of the Gribov-Zwanziger Lagrangian, see 
[U [5], within the previous set up. We start with the Yang-Mills action, as given in eq.([3]), and we add 
a set of fields {^1f,^1f), iplf^^^'^) transforming as BRST doublets. For the BRST invariant action we 
write 

So = Sym + Sgf + s J d^x (a;-a, (d^^l' + gP^^A^^^'^) - ^^^i^V^) ' (15) 



with 



sAl = -{D,cT , 

sc'' = ib" , = , (16) 
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and 



suf^ = r;, sr; = o, (17) 

sSo = 0. (18) 

For the explicit expression of Sq one gets 

So = Sym + J d\ (z6'^a^^«+c"a^(L>^c)") 

/^'^^ i „ ( Pi T-;<IC^ fabm/n „\b ,„mc , ,,2 (—ac,„ac -r-,ac, ,ac\\ (^n\ 
d X [g {d„uj^, ) } [D^c) <y2^ + ((^^ </?^ - o;^ o;^ j j . (19) 

The fields {^'Jf, ^'Jf) are a pair of complex conjugate bosonic fields. They have dimension one and zero 

ghost number. Each field has 4 (A'^^ ~ components. Similarly, the fields (ZJ^ , co'j^'^'j are anticommuting, 
having ghost number (—1, 1), respectively. 

From the BRST doublet structure of the additional fields {ip'^^ , Lp*^^ ,11]'^^ , ulf) it follows that the ac- 
tion Sq is equivalent to pure Yang-Mills theory, leading to the same ultraviolet behavior. The doublet 
structure implies in fact that there is an exact compensation in the Feynman diagrams among the bosonic 
sector (^'^, (f^) and the anticommuting one (oJ^"^, OJ^) ■ Also, it is easily verified that the integration over 
the auxiliary fields (jplf , iplf ,LJlf , cu^'^) amounts to introduce a unity in the partition function, meaning 
that the physical content of the theory is precisely the same as that of the Faddeev-Popov action. 



We now introduce a soft breaking of the BRST symmetry in a way which modifies the low energy 
sector, while keeping the original ultraviolet behavior as well as the renormalizability of the theory, as ex- 
plained in the last section. We thus introduce a quadratic term describing the coupling among the gauge 
field and the bosonic fields {^1f,^1f)- This term will affect the gluon propagator (^^(/c)A[^(— /c)) 
in the infrared. As a consequence, the correlation functions of the gauge invariant quantities, like e.g. 
(F^(a;)F^(y)), will get modified in the low energy region. For this quadratic term, one can write the 
following expression 

^7 = -I'a J d^x [r'^Al^l^ + r'^All^^) . (20) 
where 7 is a parameter with mass dimension one. With the introduction of the term the action 

S = So + S^ (21) 
does not display anymore exact BRST invariance. Instead, one has the softly broken identity. 

sS = 72A^ , (22) 

where the soft breaking term is given by 

= 1 d^gr'^ ((i^r^™) {ft ) - • (23) 

Looking at the gluon propagator {A'^(k)A!l{—k)), one finds 
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or 



safe / e kf^ki, 



{Ai(k)AU-k)) = - -I^J (25) 

with 

M^(k') = ^ . (26) 

We see therefore that, while keeping the usual ultraviolet behavior ~ l/fc^ at very high momenta, ex- 
pression (|24p turns out to be deeply modified in the infrared region A; ~ by the presence of the soft 
parameters (7, /x). The propagator (f24l) is not of the Gribov type, as it does not vanish at the origin, 
thanks to the presence of the parameter // [HIS]. Remarkably, this propagator is in good agreement with 
the most recent lattice numerical simulations [H [71 [U [9] . In particular, it exhibits violation of reflection 
positivity, a feature which is usually interpreted a signal for gluon confinement. 

Without entering into details, it is worth mentioning here that the action S of expression (j21|) has a 
profound geometrical meaning \10\ lllj . as it allows to restrict the domain of integration in the path 
integral to the so called Gribov region Q, defined as the set of gauge fields fulfilling the Landau condition 
and for which the Faddeev-Popov operator Ai"-'^ = —d^D'^ is strictly positive, namely 

n = d^A-^ = 0, M^' = -d,Df > 0} . (27) 

The restriction to this region is needed in order to properly take into account the existence of the Gribov 
copies. In particular, the parameter 7, known as the Gribov parameter, is not a free parameter of the 
theory, being determined in a self-consistent way through the gap equation [lOl [11] 



6r 

^7 



2 - , (28) 



where F stands for the IPI generating functional. This equation enables us to express it in terms of 
the coupling constant g and of the non-perturbative invariant scale Aqcd- In much the same way, the 
parameter fi has been determined by means of a variational principle, see refs.|4| |5j for the numerical 
characterization of (7, It is worth mentioning that the introduction of the BRST invariant mass term 
in expression ([19]) 

j d4x/i2(^Vr-^M'0 ■ (29) 

follows from the observation that the dimension two condensate {"^f^^lf — ^If^ff) has a non- vanishing 
value. In fact, as shown in [11[5], one finds 



/—ac,ac T-;ac, ,ac\ S^N 1 



' ' 6Att 

= 2g^Nj\ (30) 



from which one sees that (jp^^flf — lJ^oj^^ is nonzero for non-vanishing Gribov parameter 7. Let us also 
mention that the appearance of the BRST soft breaking term can be traced back to the properties 
of the Gribov region. It arises as a consequence of the fact that any infinitesimal gauge transformation 
of a field configuration belonging to Vt gives rise to a configuration which lies outside of [5] . 

A remarkable aspect of the action {Sq + S^) is its multiplicative renormalizabilit}E] [HI [HI [U [5]. 
Only two independent renormalization factors are needed to account for all ultraviolet divergences. It 



^We remind here that vacuum terms of the kind J d'^x4:{N^ — 1)7'' have to be properly taken into account when solving 
the gap equation (|28p . These terms show up in fact during the renormahzation procedure, as it follows by employing the 
method of the external sources, as outlined in [lUI 1111 [T^ HI [5]. 
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should be noticed that both the Gribov parameter 7 and the dynamical mass n do not renormalize 
independently, the corresponding renormalization factors being given by suitable products of the two in- 
dependent renormalization factors of the theory which can be chosen as being the renormalization of the 
coupling constant Zg and of the gluon field Za- The soft breaking does not invalidate the existence 
of the Slavnov-Taylor identities for the renormalized IPI generating functional T, namely 

S{T) = [A^ • r] , (31) 

where [A^ • F] is the generator of the IPI Green functions with the insertion of the breaking A^ and 

with fi^ and L"" being the external sources coupled to the nonlinear BRST transformation of the gluon 
field A'^ and of the ghost field c"". Equation (j3ip is the quantum generalization of the broken identity 
(|22p . It is worth emphasizing that the Slavnov-Taylor identities pip have in fact a powerful predictive 
character. They allow us to establish the relationships among the various Green functions of the theory 
in a way which takes into account the presence of the Gribov horizon [5j. 

Let us also mention that, due to the presence of the soft breaking, the mass parameter /u^ acquires 
the meaning of a physical parameter of the theory, even if it has been introduced through an s-exact 
term, eq. fllSp . This important feature has a simple understanding by noticing that, in the absence of the 
Gribov parameter, i.e. 7 = 0, the contributions of the two BRST doublets, (^^'^, </'^'^) and (cJ^'^, w^'^) , 
compensate each other in an exact way. In fact, setting 7 = 0, one immediately realizes that the gluon 
propagator in eq. ()24p gets independent from the parameter jj?. However, in the presence of the breaking, 
i.e. for 7 7^ 0, the compensation among the two sectors {Tp^, cp'Jf^ and (tJ"'^, w^'^) does not hold anymore, 
due to the fact that the fields (^^'^, (p'^) have been coupled to the gluon field through the quadratic term 
of eq.(|20p. As a consequence, the parameter fj? enters in a nontrivial way the gluon propagator as well 
as the correlation functions of the theory. A more formal proof of this property can be given by making 
use of the renormalized Slavnov- Taylor identities ()3ip . 



4 A model for the dynamical quark mass generation 

As second example, we present a possible mechanism for introducing a dynamical quark mass. To some 
extent, this is a way to push further the gate opened by Gribov-Zwanziger in the gluon sector, using 
the idea that all partons of a confining and asymptotically free theory should get a propagator with a 
momentum dependence as in eq. (|24p . since they all have very analogous physical status. The mechanism 
will consist in introducing BRST doublets in correspondence with the quarks, and check the possible 
BRST symmetry breaking soft terms that are allowed by power counting. 

We thus supplement the Yang-Mills action of expression ([3]) by the fermionic matter term 

5*^ = y d'^x ■07^L>^V , (33) 



with 



sS^ = . (34) 
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Further, we introduce two BRST doublets of spinor fields, (ipp^up) and {Xp,r]p) transforming as 

s(pp = ujp , sujp = , 

srjp = \p , sXp = . (35) 

The fields {ipp,up) have dimension 1/2 and ghost number (0, 1), respectively, while (Xf^Vf) have dimen- 
sion 3/2 and ghost number (0, —1). 

According to the set-up of Sect. 2, the invariant action describing the propagation of these fields is obtained 
from an exact BRST variation, namely 

S^x = s j (fx {-r]pd'^ipF + Tppd'^rjp + n?{r]pLpp - Tppijp)) 

= J d^x {—Xpd'^ipp — Tppd'^Xp — rjpd'^ujp + iJpd'^rjp 

+m'^{Xp(fp +lfpXp +fjpujp -Uprip)) . (36) 
For the soft breaking term of the BRST symmetry we get 

Sm = Ml j (fx {ippi; + i^ipp)-M2 j (fx (Xpij + i^Xp) . (37) 



From expression (|37p one realizes that, in a way similar to the case of the Gribov-Zwanziger action, 
the additional fields ifp, Xp couple linearly to the matter field ip. As a consequence, the propagator 
(^'tp{k)ip{—k)'j of the matter field will get deeply modified in the infrared by the presence of the soft mass 
parameters m, Mi, M2- In fact, one finds 

where the mass function ^(A;) is given by 

Aik) = . (39) 



fc2 + 



m 



The function A{k) is completely analogous to the function Ai'^{k'^), which has appeared in the propagator 
for the gluon in eq. pSh . Its geometrical interpretation is however lacking, since the Gribov phenomenon is 
left untouched by the quark dependance. For the quarks, the demand of a modification of the propagator 
is thus motivated by their confinement and by the chiral symmetry breaking, and it is made possible in 
quantum field theory by this general mechanism involving the introduction of soft breaking terms using 
doublets. 



Interestingly, expression (j39p provides a good fit for the dynamical quark mass in the infrared region, 
as reported by the lattice numerical simulations of the quark propagator in the Landau gauge, see for 
instance [131 [U]. 

Although being out of the main goal of the present letter, let us mention that the renormalizability 
of the spinor action (5^a + 5'*/) stems from the fact that the fields (ifp, Af) > (^f^Vf) are introduced in 
a way which enables us to write a local set of linearly broken identities, as expressed by 

= -O LOp + m ujp , (40) 

6{S^x + Sm) ra 2 lA-,-. 
= o rjp — m rjp , (41) 

OLOp 
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^ {S^x + Sm) 
S (S^A + Sm) 



ifp + m^ipp — M2V' 



(42) 



(43) 



Thanks to the fact that the right hand sides of eqs. (j40p — (j43p is Unear in the fields, it turns out that 
these equations can be converted in a set of powerful Ward identities which constrain very much the 
possible counterterms which can show up at quantum level, ensuring that the action {Sy^x + Sm) is 
stable against quantum corrections. The reader might have in fact remarked that equations ()40p — (j43p 
follow from the fact that the quadratic terms in expression (1361) contain simple space-time derivatives, 
e.g. rJpd'^ipF, and not covariant derivatives of the kind rjpD'^ipp. This term would jeopardize the linearity 
in the quantum fields of eqs. (j40p — ()43p . as quadratic and cubic terms in the fields would be present. 
These terms have to be treated as composite operators, requiring a proper renormalization procedure. In 
summary, eqs. (|40p — (j43p would loose their usefulness. Moreover, as much as in the case of a c/)^— theory, 
the presence of a term of the kind r]pD'^ipp would require the introduction of a quartic spinor interaction 
like p{r]p(pp) (Tpprjp) with its own dimensionless parameter p, implying that we would be leaving the 
space of soft parameters. Let us finally mention that a similar set of Ward identities can be derived in 
the case of the Gribov-Zwanziger action [lOl [TTl [12l SI [5]. 

5 Conclusion 

In this letter we have illustrated a proposal in order to introduce non-perturbative infrared effects leading 
to a modification of the long distance behavior of the correlation functions of the theory. 

The whole idea relies on the introduction of additional fields giving rise to BRST doublets. These 
fields do not modify the set of original observables of the theory, identified here with the cohomology 
classes of the BRST operator. Subsequently, a set of soft parameters are introduced through terms which 
give rise to a soft breaking of the BRST invariance. The simplest way to do this is by demanding that 
the doublet fields couple linearly to the original fields of the theory, so that the soft breaking is given by 
terms quadratic in the fields. As a consequence, the propagators of the original fields of the theory get 
deeply modified in the infrared. In turn, the correlation functions of the local observables of the theory 
gets modified too in the low energy region. 

We have pointed out that such kind of coupling enables us to write down suitably linearly broken Ward 
identities which can ensure that the theory remains renormalizable in the presence of the breaking. More- 
over, we have seen that the presence of this breaking allows us to give an elementary algebraic proof of 
the fact that the corresponding soft parameters are physical parameters and not unphysical ones. 

The Gribov-Zwanziger action implementing the restriction in the Feynman path integral to the Gri- 
bov region $7 has been revised within the present set up. Finally, we have proposed a model which might 
be of a certain interest in the study of the dynamical quark mass generation. The details of the renor- 
malizability of this model as well as a discussion of possible gap equations allowing to give an estimate 
of the parameters m. Mi, M2 will be the object of a more extended work |15j . 
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